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Longitudinal and transverse behavior of magnetic field of relativistic non- 
linear three-dimensional plasma wave is investigated. It is shown that the 
magnetic field of the wave is different from zero and performs higher fre- 
quency oscillations compared to the plasma electron frequency. An increase 
in the nonlinearity leads to strengthening of magnetic field. The oscillations 
of magnetic field in transverse direction arise, that caused by the phase front 
curving of nonlinear plasma wave. The numerical results well conform with 
predictions of the analytical consideration of weakly- nonlinear case. 
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The progress in the technology of ultrahigh intensity lasers and high cur- 
rent relativistic charged bunch sources permits the use of laser pulses 1 or 
charged bunches 2 for excitation of strong plasma waves. The excited plasma 
waves can be used, for example, for acceleration of charged particles and 
focusing of bunches. 2 The amplitude of longitudinal electric field E max of rel- 
ativistic plasma waves excited in cold plasma is limited by the relativistic 
wave-breaking field 3 E rel = [2(7 - l)] l/2 E WB / (3, where 7 = (1 - (3 2 )~ 1/2 is a 
relativistic factor, (3 = v p h/c is a dimensionless phase velocity of the wave, 
E WB = m e 0J pe v p h/e (E WB [V/cm\ m 0.96n,!/ 2 [cm -3 ]) is the conventional non- 
relativistic wave-breaking field, to pe = (Aim p e 2 /m e ) 1//2 is the electron plasma 
frequency, n p is the equilibrium density of plasma electrons, m e and e are 
the mass and absolute value of the electron charge. 

The linear plasma wave theory is valid when E maiX <C Ewb- If is well 
known that a three-dimensional linear plasma wave in a cold plasma and in 
the absence of external fields is potential. 4 In this case the magnetic field of 
plasma wave is zero. The magnetic field is absent also in one-dimensional 
nonlinear case due to the symmetry of the problem. However, magnetic 
field of nonlinear three-dimensional plasma wave with the relativistic phase 
velocity (/3 ~ 1) was not studied up to now and our aim is to investigate this 
problem. 
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We shall study nonlinear plasma waves (wake waves) excited in cold 
plasma by relativistic electron bunches or intense laser pulses (drivers) and 
suppose the azimuthal symmetry of the problem. In this case for non-zero 
components of plasma electrons momentum and electromagnetic field of the 
wave we have the following set of equations 5 ' 6 : 

'£-£-^ = °. < 2 > 

-f^f + ^o, (3, 

V ± H e + /3— ^ + & iV e + /5a = 0, (4) 
or 

N e = 1 — a — V L E r - (6) 

As usual, the Eqs. (1) and (2) were derived taking into account that the curl 
of the generalized momentum is zero, P 2 H - rotP = 0, or m our case 

In Eqs. (1) - (6) 7e = (l+P 2 2 +P r 2 +a 2 /2) 1/2 ,P z , r = P z , r lie and N e = n e /n p 
are respectively a relativistic factor, dimensionless components of velocity 
and dimensionless density of plasma electrons, z = (u pe /v ph )(Z — v p ht), a = 
nb(z,r)/n p , rib is density of bunch electrons, a = eE (z,r)/m e cuo, where E 
and u are the amplitude and frequency of laser pulse, Vj_ = d/dr+l/r. Also 
the following dimensionless variables have been used: the space variables are 
normalized on \ p /2ti = v p h/uj pe , where X p is the linear plasma wavelength, 
the momenta and velocities - respectively on m e c and the velocity of light 
and the strengths of electric and magnetic fields - on the nonrelativistic wave- 
breaking field E WB . In the general case the analytical consideration of the 
problem seems impossible. First we consider weakly-nonlinear case and then 
present numerical results. 
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In the weakly-nonlinear case we suppose (3 = 1 and use the generaliza- 
tion of the well known expansion 7 which was used to study one-dimensional 
nonlinear relativistic-plasma- waves 8 : 

u(r,z) = eui(r,tf) +e 2 u 2 {r,^) +e 3 u 3 {r,^) + (8) 
&f!/dz = 1 + eki(r) +e 2 k 2 (r) + 

where u stands for normalized values P ZyT , E z>r or Hg, e = E™ ax <C 1 is the 
small parameter, P 2 r , E z>r , H 9 <C 1, A p is the nonlinear plasma wavelength. 
Substituting expansion (8) into equations (l)-(5) we have: 

dP Z i 

-q^- ~ E *i = 5 M> (9) 

^ - En = % (10) 
dH ei dE ri _ 

~~d¥ + ~d¥ + n ~ 3 ' u { } 

BE ■ 

V ± H ei + -^ + P zl = S 4 „ (12) 
dHoi dE ri dE zi _ 

~d¥ ~ ~d¥ + ~d¥~ 5 ' 4 ' 1 j 

where subscript i denotes the order of approximation and Si_ 5) j are the non- 
linear functions of tti-2, and u±. In the linear case (i = 1) Si_5 ; i = 
and one can obtain well known solutions (see, e.g., Ref. 9): P z \ = — Psm*, 
P rl = (dR/dr) cos*, E zl = -i?cos*, P ri = -(dR/dr) sin*, and # ei = 0, 
where i?(r) depends on the radial distribution in the exciting source, i. e., 
on a(r) or a 2 (r). In the second approximation, from the condition of absence 
of resonance terms (proportional to sin * or cos *) one can find that k\ = 0. 
Then 



5 1>2 = (1/2)9(^ + ^1)/^, 

S 2 , 2 = (1/2)0^' + J*)/0r, 

S 3 ,2 = P r i(V ± P rl + <9P zl /«9*), 

^4,2 = P zl (V ± E rl + dE zl /dV), 

S 5 ,2 = 0, 
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and for magnetic field strength we obtain the following equation: 

[(d/dr)V ± - l]H e2 = dS 4 , 2 /dr - dS 3 , 2 /d^. (14) 
Solution of this equation is 

#w = /i(r) + / 2 (r)cos(2tt), (15) 
where fi, 2 (r) satisfies to equations 

(A ± - r " 2 - = (l/2)(d/dr)[R(A ± -l)R], 

( A± -r- 2 -l)/ 2 = (l/2)[(dR/dr)A ± R- R(d/dr)A ± R], 

here = V ±(d/dr) is the transverse part of Laplacian. In the i — th 
approximation H ei ~ cos^^). The dependence k 2 (r) can be obtained in 
the third approximation. So, the weakly-nonlinear theory predicts that the 
magnetic field of 3D nonlinear plasma wave is not zero, in opposite to the 
3D linear case and to ID nonlinear one, and oscillates at the harmonics of 
the linear plasma frequency; the nonlinear wavelength changes in the radial 
direction. 

We have solved Eqs. (1) - (6) numerically choosing the Gaussian profile 
of the driver both in longitudinal and transverse directions: 

A(z, r) = A exp[-(z - z ) 2 /a 2 z ] exp(-r 2 /a 2 ), (16) 

where A(z,r) stands for a or a 2 . Fig. 1 shows three-dimensional linear 
plasma wave excited by relativistic electron bunch. The linear numerical 
solution obtained well agreed with predictions of the linear theory. 2 As is 
seen in Figure 1(b), in the linear case the magnetic field excited in plasma 
is localized in the range occupied by the bunch and in the wake i?™ ax ^> 
jjrrmx correspond to the linear theory. 

In the nonlinear regime the behavior of plasma waves is qualitatively 
changed. In Fig. 2 one can see a nonlinear plasma wave that is excited by 
an intense laser pulse. The main difference here from the linear case is the 
change of shape and length of the wave with the radial coordinate r [see Fig. 
2(a)]. The magnetic field strength in the nonlinear plasma wave as shown in 
Figure 2(b) is different from zero and has the magnitude comparable with 
that of other components of the field. The magnitude of higher frequency 
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oscillations (as compared to the plasma frequency) performed by the mag- 
netic field along z grows in proportion to the nonlinearity. Such a behavior of 
magnetic field is a purely nonlinear effect. Indeed, in the linear case Hg = 0, 
i.e., the contribution of momentum components at the plasma frequency in 
the expression (7) is compensated. The nonlinearity of the wave implies a 
rise of higher harmonics in P z and P r . According to (7), the rise of magnetic 
field is due to these harmonics and this accounts for frequent oscillations seen 
in Figure 2(b), that conforms with the weakly- nonlinear analytical considera- 
tion presented above. On the other hand, the non-zero magnetic field means, 
according to (7), that the motion of plasma electrons in the nonlinear wave 
is turbulent (rotP ^ 0). The degree of turbulence (the measure of which is 
Hg) grows in proportion to the nonlinearity. 

It is easy to see that due to the dependence of the wavelength on r, the 
field in the radial direction grows more chaotic as the distance from the driver 
increases. In fact, the oscillations of plasma for different r are "started" 
behind the driver with nearly equal phases but different wavelengths [see 
Fig. 2(a)]. As \z\ increases, the change of phase in transverse direction (for 
fixed z) becomes more and more marked. This leads to a curving of the phase 
front and to "oscillations" in the transverse direction. 6 ' 10 The radial behavior 
of the longitudinal electric field and magnetic field of the nonlinear plasma 
wave is presented in Fig. 3. Qualitatively the radial dependence of the field 
differs from that of the linear case by the change of sign and "steepening" of 
fields along r and that is connected with the curvature of wave phase front. 

So, we have found that the magnetic field of nonlinear plasma wave is 
different from zero and performs higher frequency oscillations as compared 
to the plasma frequency. The latter qualitatively differs from the linear 
case and one-dimensional nonlinear one. The obtained numerical results well 
agreed with the predictions of weakly-nonlinear theory. 
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FIGURE CAPTIONS 

Figure 1. The linear plasma wave excited by an electron bunch with 
gaussian profile; «o = 0.1, a z = 2, a r = 0.5, 7 = 10. (a). The dimensionless 
strength of longitudinal electric field. 1 — E z (z) at the axis, r = 0;2 — r = l; 
3 — r = 2. (b). The strength of azimuthal magnetic field. 1 — r = 0.2; 2 
-r = 0.75; 3 — r = 2. 

Figure 2. The nonlinear plasma wave excited by laser pulse. The pulse 
parameters are : = 3.6, a z = 2, a r = 5, 7 = 10. (a). The longitudinal 
electric field E z . r = 0,2,4 and 5 in the order of magnitude reduction, (b). 
The magnetic field strength. 1 — r = 2;2 — r = 4. 

Figure 3. The radial behavior of field for the case given in Figure 2 for 
z = —25. 1 — longitudinal electric field E z {z = — 25, r); 2 — magnetic field 
H e (z = -25, r). 
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